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Abstract 

We propose to replace ordinary propagators in lattice operator correla- 
tions entering the determination of hadron masses with space-time smeared 
propagators. These are defined as the inverse of the quadratic operator in the 
fermion action times its hermitian conjugate. We obtain a cleaner determina- 
tion of hadron masses, comparable to the effect of using smeared operators. 
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The determination of hadron masses in lattice QCD simulations is based on 
the study of the correlation of two operators at different points with appropriate 
quantum numbers. Due to translation invariance the correlation depends only upon 
the relative distance: after integrating over the space components, it depends upon 
the time difference and, when this becomes sufficiently large, it decreases exponen- 
tially with a slope equal to the mass of lightest hadron with the specified quantum 
numbers. The expression of the correlation at a generic time can be obtained by 
inserting a complete set of states: 
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where m is the mass of the lightest intermediate state. The quality of the correlation 
depends upon the matrix element of the chosen operator between the vacuum and 
this state: it is good when the contribution of this state dominates, leading to a 
precocious asymptotic behaviour. With pointlike operators it takes several lattice 
points before the asymptotic behaviour is established, restricting the determination 
of the hadron mass to a region of large times where the signal is very small. 

Various improvements have been proposed in the literature [1], all based on 
the use of smeared sources obtained as a superposition of operators where the fun- 
damental fields are taken in a smearing region with weights which can be optimized. 
Assembling coloured fields at different points requires a parallel transport to preserve 
gauge invariance or fixing the gauge. In the first case one has some degree of arbi- 
trariness and time spent in the construction of the sum of paths for the transport, 
in the second case some noise is introduced by the gauge fixing procedure. 

In this letter we propose a new solution which makes use of pointlike operators 
but modifies the propagators which appear in the correlation after integration over 
fermions. For example, the determination of the pion mass can be obtained by 
choosing the operators 0\ and O2 equal to ^>75^> in eq. 1 which, after fermion 
integration, becomes: 

C(t) = j d 3 x(0\T${0) l5 i>{0)j>{x) l5 i>{x)]\0) 

= JvAe- 3 ^ J d 3 xTT{G(0,x) l5 G(x,0) l5 } (2) 

where G are the fermion propagator in a fixed gauge field configuration {A} and 
only the flavour non singlet correlation is retained. 

The idea is to replace the usual propagator with one which leads to the same hadron 
spectrum but which is more suppressed with increasing quark off-shellness. Excited 
states are expected to be correspondingly more suppressed. The standard propaga- 
tor is the inverse of the quadratic operator appearing in the Wilson action: 
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Figure 1: The effective mass for the pion correlation in the free quarks case. The 
data + correspond to the smeared case. 

The new propagator is taken as the inverse of the operator M ] M and can be seen 
as a smearing, also in time, of the original propagator. 

For free quarks in the continuum this just corresponds to replacing the usual eu- 
clidean fermion propagator 



m 



p 2 + m 2 



with 



p 2 + m 2 



(4) 



In this case it is easy to obtain from eq. 2 the analytic form of the correlation which 
is expected for large times to be dominated by the lowest energy state of two free 
quarks of mass m: 

C(t) oc r dE E OL \[E 2 

J m 



m 2 e 



-2Et 



(5) 



The value of the parameter a is 1 in the standard case and —1 with the smeared 
propagator, leading in the latter case to a stronger suppression in the momentum 
integral over off-shell intermediate states. In figure 1 it is shown the effective mass 
on the lattice for the two cases for free Wilson fermions. 

For the interacting case, before proceeding to a comparison of the correlations, 
we want to show that the ones obtained from smeared propagators are indeed dom- 
inated for large times by the lowest lying state. 
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The proof uses the following relation between the standard and the smeared propa- 
gator: 

(M^M)- 1 = Jlf-^jJIf 1 ^ (6) 

The expression for the correlation with meson sources and smeared propagators can 
be rewritten as: 



J d 3 x J d A yd A z M-Ky^M^yo^YM^o^M- 1 ^! 

where for pion T = 75, and reinterpreted as the following new correlation: 
J<Px Jd 4 yd 4 z{0\T[ir 12 (y)Y 23 (0)ir 34 (z)Y 41 (x)] |0) 



(7) 
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where the quark fields entering the correlation are labelled with a new flavour index; 
each new flavour has the same mass, standard propagators, but does not give rise 
to additional fermion loops and it must be considered quenched. The distinction 
between flavours is needed in the Wick expansion to produce only the contractions 
represented in eq. 7. The spectrum contains new states exactly degenerate with 
the standard ones. By summing over all different time orderings, using translation 
invariance and inserting a complete set of states one obtains the final expression for 
the continuum case in a infinite volume: 
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The main difference with respect to the usual correlations is that the parity of the 
intermediate states is not determined by the one of the operator chosen for the corre- 
lation: both parity states contribute because of the insertion of the parity commuter 
operator ip*y s ip. For large times the lowest lying state dominates: however, in the 
case where the source operators are of even parity, the coefficient of the exponential 
of the pseudoscalar state is the sum of a constant term and of one increasing linearly 
with time and leading to a slow drift of the effective mass calculated assuming a 
constant coefficient, which reaches the true value only asymptotically. This is not 
the case for parity odd operators, where such a behaviour affects the exponential 
of the heavier scalar state. In general, such a pathology occurs when there are two 
degenerate intermediate states in the particular time ordering where the parity com- 
muting operators occurr at times between those of the two source operators. 

In figures 2 and 3 we compare the it and p effective masses 1 obtained with stan- 
dard and smeared propagators on a 16 3 X 32 lattice at (3 = 5.7 for quenched Wilson 
fermions at K = 0.161: in the latter case the plateau appears at smaller values of the 
time distance. The points referring to scalar and pseudovector sources are obtained 
from a naive analysis ignoring the pathology of the term linear with time. A proper 
fit which allows for this behaviour leads to values of the mass compatible with those 
of the standard and smeared propagators with pseudoscalar and vector sources. 

For baryon states, the technique which allows to show that the correlation is 
indeed dominated by the lowest lying state of any parity is similar to the meson 
case. The correlation of the two proton operators of the form: 

f d 3 x{0\T[y 6 (x)V 6 '(0)}\0) 
with V s (x) = [ U T{x){C l5 )d b {x)] u s c {x) e abc , 

$ = $^7o and C the charge conjugation matrix (10) 

can be transformed into: 

Jd*xJd i yd i zd i w(0\T x^i/K^K^) (* 2 (0) 75 f ] |0) 

where = [u T a .(x)(C l5 )d bii (x)} u s c -(x) e abc 

and tt£ = fa 5 fj (11) 

using the relation in eq. 6 between the space-time smeared propagators and the 
standard ones. By inserting a complete set of states and performing the integrals over 
times as in the meson correlation case one gets a sum of contributions from baryons 
of both parities. An interesting property of such a correlation is the symmetry with 
respect to the center of the lattice, which allows to further average statistically the 
signal. This can be shown by using the representation of the correlation in terms of 

1 The effective mass is obtained from an analysis of the correlation which assumes only one 
intermediate state at any time. 
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propagators with their periodicity properties and by considering the contributions 
coming from the correlations forward and backward in time. 

In the case of baryons, it is not possible to eliminate by a suitable choice of 
the source the contribution where there are two degenerate intermediate states with 
the lowest mass because of the insertion of three parity commuting operators. The 
determination of the effective mass is then biased, at moderate times, by the presence 
of a coefficient linear in time, although in this case one gets a smaller percentage 
correction with respect to the pion case. 

A solution to this problem is the choice of a single smeared propagator and two 
standard ones. The symmetry of the correlation is lost in this case, but it can be 
shown that, for the best case where non relativistic sources are used, the signal from 
the correlation backward in time is very small and leaves a sizeable "fiducial" region 
for the determination of the effective mass plateau. 

The final results for baryons are shown in figure 4. In the first plot we show 
the effective mass for the following non relativistic proton operator: 

*i(x) = [u T a (x)(P ± C l5 )d b (x)} u s c (x) e abc with P ± = 1(1 ± 7o ) (12) 

with standard propagators on the up quarks and a smeared one on the down quark 
compared with the cases with three standard or three smeared propagators. In 
figure 5 we show the absolute value of the real part of the correlations for the same 
three cases. The points where the correlation changes sign are affected by large 
statistical fluctuations and are not shown in the figure. 

The use of space-time smeared propagators has been shown to improve signifi- 
cantly the quality of the correlations from which hadron masses are determined. The 
extra computing effort is rather limited and needed anyway if the conjugate gradient 
algorithm is used in the inversion: standard propagator correlations are obtained 
by smearing the quark fields in the operator creating the initial state with the local 
operator . The technique presented in this letter can be further improved with 
the standard operator smearing. 
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Figure 2: The pion effective mass on a 16 3 X 32 lattice at (3 = 5.7 for three differ- 
ent cases: (□) standard propagators, (o) smeared propagators with a pseudoscalar 
source operator, (o) smeared propagators with a scalar source operator. 
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Figure 3: The p effective mass for three different cases: (□) standard propagators, 
(o) smeared propagators with a vector source operator, (o) smeared propagators 
with a pseudovector source operator. 
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Figure 4: The proton effective mass for the cases: (o) three standard propagators, 
(o) one smeared and two standard propagators, (□) three smeared propagators. 



le+10 



|ReCp| l 



le-10 




10 



15 



20 



25 



30 



Figure 5: The absolute value of the real part of the correlation for the three cases 
defined in the previous figure. Points where the signal changes sign and is affected 
by large statistical error are omitted. 



8 



